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1. Introduction 
This final report summarizes the work performed under Contract 
No. F19628-78-C-0206 which the Electronic Systems Division of the 


Air Force Systems Command granted to the EMtec Engineering, Inc. 


completed in August, 1979. 

The major objectives of this R §& D study were to perform 
theoretical and numerical analysis of multimode fiber optic components. 
Specifically, the following multimode components were studied: (a) 


Dual Fiber Couplers, (b) Fiber Tapers, (c) Fiber Horns and Branching 


i 
; 
| 
Los Angeles, California. The work was begun in August, 1978 and | 


Waveguides, and (d) Mode Converters. In the following we shall 
first present some background information on this subject and then 
our mathematical approach will be delineated. In Section III 
selected results of our study are then summarized. Finally, the 
concluding remarks and recommendations for future work will be 
given in Section IV. Listings of the important computer programs 


that we developed are also included in the Appendix. 


II. Background and The Mathematical Approach 


There are currently strong trends in system design toward 
microminiaturization digital processing and system level integra- 
tion in order to achieve smaller size, weight, and power consump- 
tion, along with lower cost and improved reliability. These trends 
naturally point to data bus multiplexing, i.e., the interconnection 
of a number of spatially distributed terminals via fiber optic 
waveguides cables. The key component for any fiber data bus system 


is the fiber coupler. Since multimode singlestrand fiber is used as 


a communication link in a data bus system, multimode fiber couplers 
must be designed and used, 

To design any coupler properly for a multimode single fiber 
data bus system, detailed analysis of waveguide propagation must 


be carried out. Existing techniques based on the finite elements 


method, ! the coupled mode theory, or the geometrical optics method” 


are usually inadequate. (Detailed discussions have been included in 
a review paper which would appear shortly.) Although the finite 
elements method is a very powerful approach when dealing with single- 
-mode fibers or couplers, it is very inefficient and costly (in terms 
of computer time) to obtain any results for multimode structures. 
Similarly, coupled mode theory has been used quite successfully in 
predicting the coupling efficiencies of single-mode structures, but 
cannot be used for multimode structures. Although the geometrical 
optics method using the ray-tracing technique may yield zero-order 
results for multimode structures, it is too crude to predict the 
wave behavior of light signals in couplers. Recently, Arnard> 
developed a technique based on the Cook adiabatic coupler principle® 
to treat the multimode coupler problem. He claims that "it may be 
used to couple two optical fibers because the dimensions are not 
critical; only slowness is required". Our technique, given below, 
will provide more accurate results because it does not involve making 
a priori assumptions on the relevance of the dimensions of the 
structure and the slowness of the coupling. Furthermore, this 
technique may be used to study mode conversion effects due to the 
nonuniform distribution of dielectric material withint the guided 
structure. 


It is noted that analysis based on the scalar wave approxima- 
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tion has been very successful in predicting the properties of many 
optical devices whose characteristic dimensions are on the order of 
many wavelengths. The scalar wave approach is therefore preferred, 
We have developed an efficient numerical technique based on the 
scalar wave equations to solve problems dealing with multimode 
couplers. 

It is well known that solving the exact electromagnetic 
equations for a spatially inhomogeneous medium is a formidable 
ech iem. Multimode couplers and mode converters can be approxi- 
mated by structures with a spatially nonuniform dielectric mediun. 
Fortunately, there are some approximations that can be made to 
simplify this task. First, the light wavelengths of interest are 
much shorter than the inhomogeneity scale length. This enables 
us to neglect polarisation effects. Hence, the optical field can 
be derived from a scalar u(x,z) that satisfied the reduced wave 


equation’ 
(v7 + un? (x,2)} wlx,z) = 0, | (1) 


where k is the wavenumber 27/A, X is the laser wavelength, and 
n(x,z) is the spatially inhomogeneous refractive index of the 
medium. Second, if we write u as the product of a factor ert, 
that accounts for the rapid change in the phase of u along the 


direction cf propagation and a complex amplitude A(x,z), a further 


simplification of the calculational problem results: 


2 
' ) z t by at . - _ 9 A(x, 2) 
[i2kn, $s + Vy +k{n (x, z) - Ny A(x, z) ae eh are, } (2) 


2 + a*/ay? and ny isa 


where ve is the transverse Laplacian 92 7ax 


given constant which represents the refractive index of some uni- 
form medium. At laser wavelengths, the complex amplitude A(x) 
varies much more rapidly transverse to the direction of propagation 
than it does along the divection of propagation. This enables 
us to make the paraxial approximation and neglect the term on 
the right side of Eq. (2) (in the Russian literature this is called 
the parabolic approximation). So, the complex amplitude now satis- 
fies 

2 ne 


+ 2 + 2. +e (n°(x, z) s 2) A(x,z) = 0... (3) 
ox ay ab. 


Z 
In addition to Eq. (3), the complex amplitude satisfies an initial 
condition on the fiber end; at z = 0, 
A(x,0) = u(x,0) 
and the boundary condition is 


A(t, et (5) 


If a truncated guassian beam is focused on one end of the optical 


guide, then 
2, 2 age 
u(x, y, ©) Us exp(-r /w ) for0 £r<b 


és (6) 
= Q for fob, 


where x? = x? + y’, w is the spot size of the beam, and b is the 


radius of the truncated beam at z = 0. 


We have solved Eq. (3) numerically in the following manner. 


Let us write A(x,z) in the form 
A(x,z) = exp[T(x,z)] v(x,z) , (7) 


where ['(x,z) is a phase function associated with the medium inhomo- 


geneities 
A 


I(x, z) = 2n, [n?cx, y,z') - no | dz! a (8) 


2 
oO 


The modified complex amplitude v(x,z) then satisfies the equation 


r 


i2kn, 2. v(x,z) te” vifet v(x, z)] = 0 (9) 


dz 


Although Eq. (9) does not look any easier to solve than Eq. (3), it 
is easier to solve numerically because, for sufficiently small 
increments in the z direction and an Bhnoneeaiely chosen lower 
limit in the integral in Eq. (8), the value of v(x,y,z + Az) can 


be obtained to a good approximation by solving the simpler equation 


: a 2 
[izin, SS oP ve v(x, io mit 6 aes (10) 


with the initial condition 
v(x,y,0) = u(x,y,0) (11) 


Physically, these equations approximate the propagation in the 
inhomogeneous medium by a two-step process at each z increment. 
First, we propagate the field u(x,z) at z to z + Az assuming that 
the intervening space is homogeneous. The effect of the inhomo- 
geneities between z and-z + Az is then accounted for by multiply- 
ing this solution by the phase factor exp!. 

In this research we hava solved Eq. (10) by the fast Fourier 
transform technique. Replacing the Laplacian by its finite 
difference equivalent but still retaining the z derviative, the 


solution of Eq. (10) can be expressed in the form 


N-1 
v(m,n,z) = p2ee V(m!',n',z) exp AST (mm! + nn')] ; 


m!',n'=0 
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where x = mAx, y = nAx, and V (m', n',z) satisfy 


iZkn 2. , Amn) V(m',n',z) = 0 
igh] {Ax} 


with the initial conditions 


N-1 
toy = = y 
¥im" in"; z.) = 3 : v(m,n, z) exp [T(m, n, z-)] 
m = 


exp]- 32% (m'm + n'n) | S 


The function f(m',n') is determined by the difference approximation 
used to represent the Laplacian in Eq. (10). For example, if Vav 


is approximated by the simple central difference expression 


ve — [v(m + 1,n,z) - 2v(m,n, z) 
az - 2 
(Ax) 


+ v(m - 1, n,2z) ren ee + 15%) 


= 2v(m,n,z) + v(m,n - 1, z)] 


then f(m',n') is . 


: ; 
£{m',n') = -4[ sin Ce) + sin? (22)| sit (16) 
Note that the series in Eq. (12) is simply the discrete 
Fourier transform of the function V(m',n',z), and thus can be 
evaluated numerically for a given V(m',n',z) by a fast Fourier 
transform algorithm. Furthermore, the functien V(m',n',z) is 


readily determined from Eq. (13) as 


feifon', nv) |. ) 
V(m',n',z) = V(m',n'z.) exp (Ax) “no 


at + cdldhaiasiabtl 


where v(m',n",z;) is given by the series in Eq. (14), which can 


also be evaluated by a fast Fourier transform algorithm. To 


summarize, we will step from z to z + Az as follows: (1) take the 
inverse discrete Fourier transform of u(m,n,z) = [exp I(m,n,z)] . 
v(m,n,z) by means of an inverse fast Fourier transform algorithm; 

(2) multiply the result by exp [-if(m" yn") A2/2k(Ax) 7n J] and take 

the discrete Fourier transform with a fast Fourier transform 
algorithm; and (3) multiply the result by exp[I'(m,n,z + Az)] to yield 
u(u,n,z + Az). This process is repeated until we have reached the 
desired z plane. 

Using this algorithm we have been successful in obtaining many 
meaningful results for various multimode structures. These results 
are summarized in the following section. 

III Multimode Fiber Components 

Before we proceed with the presentation of the numerical results, 
it should be recalled that we are dealing completely with total field 
quantities and not with the modes. In other words, we are interested 
in how the total field evolves as it propagates down the guiding 
structure; we are not interested in how each mode propagates. Never- 
theless, it is recognized that the total field may be decomposed into 
a set of orthonormal guided modes. For example, an incident Gaussian 
beam with a given beamwidth w may excite many modes in a parabolic- 


index-profile fiber when & >1 or when a<l1 with 
a = (2a) /[nnw? (28) 2/7] 


or may excite only one mode when @ = 1. Only the a = 1 single mode 
will propagate down the parabolic-index-profile guide without experi- 


encing the focusing and defocusing effects. For a # 1 cases, the 


Ait eee tt HOE eT Soe gedit nk 
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input Gaussian beam will experience focusing and defocusing effects. 
Another way to interpret the above phenomenon is that multimodes 
with different propagation constants are excited by the input beam 
when a # 1, while only single mode is excited when a = 1. The ‘ 
stronger are the focusing-defocusing effects, the higher is the 
content of different modes. 

(a) Dual Fiber Couplers 

The geometry of a multimode inhomogeneous fiber coupler is 


shown in Fig. 1. Two graded-index fibers with index variation given 


ee woe: 


by 


r 
Z Tes Oe; 
n(ry_2) ~ Bs ae : , 
Lz 
where r, 5 are the radial coordinates of 1 and 2 fibers, respec- 
> 
tively, and n,, 6, a are all known constants, are fused together 


as shown. The separation distance between the centers of the 


fibers is d. A Gaussian beam representable by 


2 
u(x,y) = Uy eXP | (« + ‘| . yf ? 


where u(x,y) is the scalar wave function of the beam, and uy» W 
are given constants, is incident on one of the fibers. We wish 
to learn how this beam evolves as it propagates down the coupled 
structure. In other words, the coupling distance and the beam 
shape will be obtained. 

Results of our investigation on the multimode dual fiber 
couplers are shown in Figs. 2-5. It is seen that if two identical 


parabolic-index fibers were placed side by side with each other 
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and if the beamwidth were so chosen that a = 1 is obtained, one would 
expect the guided power to interchange between the two fibers in a 
periodic manner with a distinct single coupling length. On the 

other hand, if a departs from unity, many modes are generated and 
many beat coupling lengths occur, so the guided power will no longer 
interchange among the guides in a simple manner. For the a<<l or 
a>>l case, even larger numbers of modes are excited. Depending upon 


how the incident energy is distributed among the excited modes the 


back and forth (periodic) power exchange phenomenon could still 


prevail if the major amount of incident energy is distributed in 
several low-order modes. Only when the incident energy is distributed 
evenly among all modes and only in the limit of infinite number of 
modes will the input power be split in a 50-50 manner between the 
two fibers. 
Sample program listing and printout for the evolution of field 
intensities in the two unequal fiber coupler are given in Appendix A. 
(b) Fiber Tapers 
Fiber tapers have been commonly used to inter-connect two 
optical waveguides of different cross-sectional sizes or shapes. 
The effects of the transition on the guided fields are shown in 
Figs. 6-10. It can be seen that very profound effects on the 
guided fields are found. Extensive mode conversion may be present 
even for gently varying tapers. Sample printout for the evolution 
of field intensiters in the presence of the fiber tapers is given 
in Appendix B. 
(c) Fiber Horns and Branching Waveguides. 


Numerical computations to obtain the field behavior in a 
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multimode fiber horn structure or in a branching fiber guide have 
also been carried out using the similar technique outlined earlier. 
The cross-sectional index profiles are specified for each step in 
the direction of wave propagation. Shown in Figs. 10-11 are plots 
of beam waist as a function of the longitudinal distance for various 
different incident beam sizes. 

(d) Mode Converters. 


There exist many occasions when the presence of high order 


modes rather than_low order modes is preferred in a multimode 


fiber. It is known that nonuniform distribution of dielectric 

materials with different indices of refraction within a fiber may 
cause mode conversion. Taking advantage of this phenomenon, one 
may design a low-order-modes to high-order-modes converter using 
nonuniformly distributed dielectrics. For example, one such con- 


verter may be the following structure: 


The cross-sectional index profile may have the following form: 
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It is anticipated that the low order modes will be affected more 


by the discontinnuities. Thus energy in the low order modes may 
be converted into that of the high order modes as expected. Other 
means of generating high order modes using variations of the above 
scheme may also be envisioned. 

To evaluate the efficiency of this type of mode converter, 
and to learn the intensity distribution of the resultant field, 
the evolution of an incident beam through this channel should be 
obtained. We have successfully developed a program which is 
capable of providing such information. Illustrative example is 
given in Appendix C. It can be seen that ring-shaped higher order 
modes can easily be generated by the introduction of a slight 
dip in the index profile as proposed above. 
IV. Conclusions and Recommendations for Future Work 

We have successfully developed a computational technique, 
based on the scalar wave - FFT method, to treat problems dealing 
with various multimode fiber components, such as, fiber couplers, 
fiber tapers or horns, fiber branches or mode converters. The 
only significant restriction to keep in mind is that index variations 
of the structure under consideration must be gentle®. This consid- 
eration is normally satisfied in most practical situations. The 
basic objective of this R & D study is, thereby accomplished. Listing 
of the computer program is given in the Appendix. By simply 
specifying the index profile at each longitudinal z-step and knowing 
the initial beam shape, one may generate the propagation characteris- 


tics of the beam as it propagates down the structure. 


Using this program, we were able to trace the evoluticn of 


a given beam as it propagates down a given multimode graded index 
fiber structure. If the multimode structure were a multimode 
coupler made with two or more paralled graded index fibers in close 
proximity of each other, we can predict the coupling distances as 
well as the power distribution of the beam in such coupler; if 

the structure were a tapered or a flared multimode fiber, we were 
able to learn whether a given tapered structure or a given horn 
structure could still confine or guide a given beam; if the struc- 
ture were a mode converter made with multimode fiber containing 

an on-axis dip in its index profile, we can learn the effectiveness 
of such a dip in converting lower order modes to higher order modes. 
To achieve such capabilities,’ an extensive amount of computer 
software was developed. 

In the course of these studies, a number of important new 
research topics were generated. For example, it is known that the 
FFT scalar wave approach requires that the index profile be gently 
varying. But the quantitative definition of "gently" must still be 
specified. (We wish to push our ability to deal with almost step- 
index profile fibers.) Initial indication (through our computer 
results) seem to show that even with rather steep index gradient, 
adequately accurate results were still obtainable. This problem 
also leads to the need for us to study the possible usage of 
adaptive coordinates to improve the beam resolution as it propagates 
down the multimode structure. Another problem which is directly 
related to the case of induced beam radiation caused by variations 


in the guiding structure is also of great interest. Perhaps the 
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introduction of a lossy surface for the outer boundary of our 
computer mesh may provide a means of estimating the radiation loss. 
Associated with beam radiation problem in the case of beam reflec- 
tion caused by variations in the guiding structure, the reflection 
coefficient may be obtained by correctly summing over incremental 
reflection coefficients (i.e., taking into account the shibes of 
the reflected waves) which occur whenever wave energy propagates from 
a region with a certain index value to another region with a different 
index value. 

With the help of our newly developed program, the principal 
thrust of the future work should be to study realistic multimode 
fiber structures, such as data-bus type multi-channel couplers, 
tapered or flared transition joints, multimode fibers with non- 


ideal parabolic index profiles and to obtain data for these realis- 


tic situations. 
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Le Ae TOT THEA ARS T DIMENSION UF MATHER AMATGS tes Tet iy, eo Peers 
Cc JOIM CTHE SECOWS DIMENSTUN OF WATRIK AMAT) Sri 
C IMiN CTrReE BEGINNING I-CQURDINATc FUR The PLOP) 
c [MAX (TRE ENDING I-CUORDINATE FUR THE PLST) 
32, eae IDEL (CUSeD Tis THe SeNSt ...Letmlvy LMAX, LSEL) PARIS E Ae aes OF > 
C JMiN, IMAA, JOEL CSEE IMiN, LAX, anv LTDEL LeoCrIFT Lous) “a 
C AMIN (The MINIMUM VALUE USED FUR SCALING--Sec DELON 
Cc AMAX (THE MAXIMUM VALUE USED FUS SCALTING=--stt vELON 
— € OL TLE (CHARACTER veCIGR OR HOLLERITH SiFiuG FOR TITLE) : 
> NCHAR CNO. OF CHARACTERS {nw Tlite vECTuUrR; Appa ey 
te SUB8PRUGRAM REFERENCES--(NONE) 
—€-COMMENTS==~ ae) 
C 1. CREATED 7/2/74 BY C. SOUTH 
e 2. NLEVEL CAN GE NU LARGER THAN 10. IF =5, CHARACTERS 
E C "ARE NOT OVERPRINTED. IF BETWEEN 0 Abs 9, CrARALTERS wil 
C UVERPRINTED NO MORE THAN UNCE. IF =1e, Tht bessest Chanact 
C WILL BE OVERPRINTCO TwICeE, 

C 3. SCALING IS LINEAR BETWEEN AMIN AND AMSx. IF AM EN=aiaax, 
SCALING WILL pe AUTOMATICALLY SET TO =ULL_RANGE Fur Tre MaleIA. | 
CT MaCHINE CUMPATIOILITY--Chanot VARLABL= sCPe IN DATA StATEMent 
c TO INDICATE NuMSER OF CHARACTERS PER aGRD FuR TRE MACHEN. 

C S. PLOT BORDER=-DEFAULT 1S A BURDER OF A5TeERiSKS. Tu ReSET THE 
Lo ee eee CA ARMC TER, CHANGE THE VARTAGLE {80RD~ fa [Ae VATA STATEMewt 
C CALGORISHM wILL CORRECTLY HaNuLE THE SLANn GORDER CaSc). 
C 6. CHANGED 7/51/74 TO HANDLE CASe UF 5S CaARS/ HURU. 
C 7. ChAWGeD 8/12/74 TU CHANGE SIXTH CHARS CTe Tu A STimMPLE ‘x! 
a CC 8B. CHANGED 8/2//74_ TO DEFINE Crar._S-7 45 eye ios Sib), PLUS 
Cc 9. CHANGED 9/i2/77 TU ALTER PAGE HANDLING. TRIS KkUUTING wILL Nua 
C TRY TO FIT SOTH THe SCALING TwFURMATISN AuwD THe GeeySueAte PLOT 
is ON TRE SamE PaGE, wiTH SCALING INFORM-TIO. scree Inc Plu. 
EG -NLMAX REPRESEWTS THe MAXITMUN RONS FOX Tre fi.Pul HATRLA St LLt 
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c ALLUN NG Tre SCALING INFURMATICN TO APPEAR UN TRE SAME PALE. ___ Grrl Si, 

Ta aL a Na id lich aaa a GN a ws I ie GRE tSe 

Cc siden (Wis GREYS: 

C SCALE, IF NECESSARY $ Girt. 

C SAM: og! SIE ee me Me WES ENS R eR eM eS. 
ae AMxX=AMAX Grftse 
[Fite= TAOS (LF ILEX) Gnev3e 
AMNNZAMIN GRE e3e | 
: IF (arAa.GI.AMINn) GO TO 100 _ Y oi BSNL a RTS eee 
ae AMN=L.EeU GREIS.. 
i BNaz-1.E20 | ; GreaS_ 
vU SU J=IMIN, IMAX, IDEL Gres. 

DO SU tegMiN,IMAReODEL 9 oe NERERS S Saere GreEYS:: 

+ pas Shins ab INT CANN, ABATET, JT) Gres.) 
Su “MASAMAXICAMK,AMATCI,J)) - GekETSE) 
Luo eon TiNue GrRets.. 
“hdres y Perinat ae Ge dL al a Se ay pallies GrEais. 
a ae PRINT HEADER Gre tse 
& GrErSc. 
WHORDS=(NCHAR@=1L)/NCPAt1 u GReEYst. 
TE QNCP ob. 6) wRETo CIF ILE, 1040) CTITLECI), Loi eWaGR03) GRETSE 
ee LFQNCPa.cO.d) write CIF ILE, 1041) CTLTLE CT), Lal, wis) Grersu. 
LF(NCP 7.10) WRITECIFILE, 1042) CTLIFLECL), L=1,%.0%035) GREYSL 
TFC(NCPa.EG.S) wRITECIFILE, 1043) CTITLECI),-1=1,iinRdS) Grers. 
>. +. JPECAMA. EG.AMN) RETURN REPT Me st CEOS Sor) Pee IRN 2h cd Mart TS, RL ger GREYS 
oy rk : “GREYS. 
c PREPAKE A LINE FUR PRINTING GRerSc 

Cc GREYS. 
__ LINE CL, 1) =1BUROR RSs toh 6 hes Sis Salton REE Ds GREITS 

eo a ge gaa gi ea GREYSS 
: LINE CJLOCt2, 1) =L150ROR GREIS. 
: LiWwetl,2)=LBLANK GREYS. 
eee Ae I NE ( 1 23 J= IBLANK rae Re iS. See Re eins a et Sek ON ee ee eR Br GreEYSe. 
Eray LINE (TILOC +2, 2) =TSCANK GRETSL, 
Lfwe (JLOC 2,35) =ITSLANK GREYS 
LASTLISJLOC+2 J GreySc 

_ HLASTL=EL4ASTLI-1 Ait CIT Sb WEE! ARAN 5d Oe he canker a hg GREYSC 

r WRITE CIFICE, 1045)~ lie T=1,CASTLI) GRETSL 
UU 400 LT=lMIW, IMAX, IDE . GREYSe 

p@ 250 Il=1,3 GabiSe 

_ 60 250 J=2,NLASTL_ A ADS aaa. ds ve ae: e __. GRE TSE 

—~ —“250— LIRE CI, LIJSTSCANK “oT a GREW st 
LEVv=l GREYSL. 

Jme : GRETSC. 

=4 vO 300 JJ=JMIN,JMAX,JDEL Ride ae See ae ia tehe St ee | ee oe GRETSC. 
i aa J=JS+1 ~~ GREYSC 
DEtAMAT (Ip JJ) “AMN)/ (AMX= AMN) *FLOAT(NLEVEL) +1. GRETSC 
L=MAX0(1,L) GREYSC 

3 _ _L=EMINOCNLEVEL,/L) eae we. Deletes hase be hN eae Ne GREYSC 
pera LEV=iAXOTLEV>LEVEUC(LID) GRETSC 
LEVNOW=LEVEL(L) , GrE¥se 

DG 300 K=1,LEViNOw Gre rst 

LINE CJ, KJSICHARS(L,+ DP Se OM ie TL REN se IE Se SM) GRE ESS 
~“5U0 ~~ CONTINUE GREYEC 

c GrRt Vou 

C FINO LAST PRINT PUSLTION .GREYSC 

a ene ws ton oe a RE TS 
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bi-é GNEYSC 76/176 UPT=2 ne eRe eee. oS Ba eget rime 
eo DE GID: RISE CRY GREYS 
vO 350 K=1,LCASTUL la is a Na moe aE | ale 
AK=LASTLinK et] GRkECSC 
EFCLINECAK,KL) NEJLTSLANK) GO TU $75 Grr yS 
Sere LONTINUE > BIO : ; GkKEYSE 
375° COmTiNue a 
a Ue -GTLO.AWDSKLIEQ LEV) nL TECLFLLEs 1050) (Line Cli-h_), GRxEYSC 
=I, Gaby St 
~-——,AF UTE ILES. el Fo Op SND. KL EU. 1) BRITECIEILE 1090) (LINECAT KL), | Bey 
» =l, GrEtr 
ater ee eGT.O.ANDSKL.NELLEV) WRITECIF (LE, L000) CLLNECiL- AL), Gey 
=1,KK) GixE 
AF CIFILEA.LT. 0.460 KLE. L) WRITECIFILE, 1000)  CLIMEC II sKL). GREY 
. LLIL=1,KK) Geter 
409 CONTINUE BxEY 
c Grey 
a a PPE CER AST Re, TCE Re ee ee hee gene ~CxEL St 
yXE XS 
. WRITECIFLLE, 1070) CISORDR, L=1,LASTLI) peer as 
RETSL 
ee ST CE ee i Oe a YL ee, ee ee ane NES," Te BREA SS 
L SKE TOL 
[F(NCPW.EQ.6) wRITECIFILE,1000) (TITLECI), L=l,wuRrcs) GRY 
TFCNCPw.eG.4) aRITECIFILE,OUL) CTLITLE CI), L=lewwuRros) GREYSL 
- {FinCPwe6C.10)_ wRITECIrFILE, 1002) CTITLeEC{),[=teheORVS) GRELS 
LFCNCPa.cOL5) ARITC CIF ILE, 1003) CTiTLECL), Lal, NiiuRvS) Gre 
wRITECIF ILE, 1010) fey Grey 
pO 200 Tot AL eveL eles 
y ee oe Sere to edie ekg se NG, SOE IATA FW Lies 
s AMIN=FLOAT CL -1) /FUCUAT CNCEVEL) * CAMX=AMN) tAMA Grey 
XMAXSXMINtDELTA GREYSCL 
ee een ey peter 
6 ot ee DUE De es ee aGREN Sle 
"gig > sa lear NOSJCEGLCEV) WRITECIFILE, 1020) TChARSCE, J), XIN, HET SE 
rR reySe 
IFCIFILEA.LT.O.AND.J.EG.1) WRITECIFILE, 1020) ICHARSCL+J)/XMIn, GReSC | 
bk RAK so, © Brak Roe er eS BEY Seo 
JECIFICES.GT.0.% AND. J. NEVLEV) WRITE CIFILE, 1030) “ICHARS(E, J) GRE 1ST 
IFCIFILEALT.OCANDCJCNESL) ARITECTFILE, 1630) Icns RSCied) GREoL 
175 COnTisve i i Greyse 
Pee SEOs TE PRT ed ee ar See ie OSU ne Ricans i op ites Le eT Be 
RETURN GReEySC 
1000 FORMAT(1HO,20A0) Grerse 
is Screg eA ed BREYSL 
he ah GRAAT Cine, L2eAl0) st ie ee Be a TE sae NO IPR eer eR eee SENS AR EOIN CAS 
eos ronmarti ao p2da5 ¥ > 5 GREYSL 
14910 FORMATCIin0, 32HGREY=SCALE CHARACTERS AND RANGES/IX) GRESCE 
LG2U FURMAT(5%,A1,5X%,E15~26/5X,E15.6) GRETSC 
BE cc hl Opi at te ge erate a lle cc pe 
1039 FORMATCIAI) GREYSe 
1940 FORMATCIH1I/1X,20A6) GRETYSC 
1241 FORMAT CIn1/1%,20A4) GREYSC 
CPGG2 PORMATUINISI As 12R1 02 oo : 2k ES 
1045 FORMATCIAL/IA, 2445) GREYS 
L045 FORMAT(1A0,1352a1) GREYSL 
$yuS0 FORMATCIA, 135201) GRE Sc 
1060 FUx%AT (int,  .152A1) oan RE YL 
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MISSION. 
of 
Rome Air Development Center 


RADC plans and executes research, development, test and 
Selected acquisition programs in support of Command, Control 
Communications and Intekligence (C51) activities. Technical 

and engineering support within areas of technical competence 
4s provided to ESD Program Offices (POs) and other ESD 
elements. The principal MLSSLON ALeas are 
communications, electromagnetic guidance and control, sur- 
veillance of ground and gekearabe objects, intelligence data 
collection and handling, information system technology, 
donospheric propagation, solid state sciences, microwave 
physics and electronic reliability, maintainability and 
compatibility. 
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